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Comments on ”Dimerization Induced by the RKKY Interaction” 


Recently, Xavier et al. claimed the existence of an in- 
sulating spin dimer State in the one-dimensional Kondo 
lattice model at quarter-hlling amidst the paramagnetic 
metallic phase [1]. Such dimer state is characterized 
by q = tt oscillations of the nearest-neighbor localized 
spin correlation, (Di) = (Si -Sj+i), under open boundary 
conditions(OBC) in their density matrix renormalization 
group(DMRG) analysis [2], where they increased the sys- 
tem size up to L= 120 keeping m = 800 States per block. 
In this comment we show that such oscillations of (Di) 
do not persist in the bulk limit and that the correlation 
function, ( DiDj ), decays algebraically at long distances. 
We first present in Fig. 1 the value of dimer-dimer cor- 
relation function, \(OiOj)\, Oi = Di — Di + \ obtained by 
the DMRG under several boundary conditions(BC). As 
for the OB C it clecreases first, and then increases near 
the boundary of the system (i.e. \i— j\ ~L) as expected 
from the report of Xavier et al.. However, when we add 
an extra site to the system and put the potentials, /i8 
and /i( 1 — 5), at the end sites, it decays signihcantly fast 
towards zero. Here, we control the potential /r to keep 
the density of the electrons at quarter-filling. We find 
that with increasing system size, these two results be- 
come close together and seem to asymptotically approach 
the line, \(OiOj)\ oc \i— j| -0 ' 5 . This power law decay is 
conhrmed by the results obtained under the anti-periodic 
boundary conditions(APBC) as shown in the hgure. 



FIG. 1. \(OiOj)\ as a function of i — j\ at J = 1 under 
OBC(* = 36, 44 for L = 48 - 49, 56 - 57, respectively) and 
APBC(L = 32). The dotted line proportional to \i — j\~ 0 ' 5 
is a guide for the eye. The inset shows the results at J = 0.5. 

We next discuss the effect of BC in detail. Figure 2(a) 
shows \(DiDj)\ when the BC is modihed from <5 = 0 to 0.5. 
The correlation function varies smoothly from the one at 
(5=0, where (5=0 reproduces the result of the usual OBC. 
The corresponding (Dj) given in the inset of Fig. 2(a) 
shows that the large oscillation of (Dj) characteristic of 
the dimer state disappears at 5=0.5. It should be noted 
that the energy density away from the boundaries is the 
same between 5=0 and 0.5 within the truncation error. 
These results indicate that the spin structure is sensitive 
to the BC when the translational invariance(TI) is vio- 


lated. Actually, (Di) under APBC and OBC at 5=0.5 are 
almost «-independent, i.e. the system has TI. Therefore 
\(DiDj)\ at 5 = 0.5 is almost the same as that obtained 
under the APBC as shown in the inset of Fig. 2(b), where 
we have removed small Adependence in (Di) by taking 
an average over the system. These results indicate that 
the dimer structure is induced by the BC which breaks 
the TI, and it is not an intrinsic bulk property. 



set shows (Dj) under OBC and APBC. (b) Charge gap, A c , as 
a function of IT 1 . (The inset shows \(DjDj + i)\ of APBC and 
the system-average of \(DjDj + i)\ under OBC with 5=0.5.) 

The proper examination of the present issue requires 
severe numerical accuracy. We calculated the system size 
dependence up to L = 240 taking a maximum of to = 1100 
to confirm the absence of a dimer long range order. In 
faet, more than to = 600 States are necessity even at L= 32 
for J= 1, and we found that one can no longer get enough 
reliable results for L>80 at J = 0.5 even with to = 1100. 
We also calculate the charge gap, A C (L), at J = 1 where 
the numerical accuracy is well controlled by to. As shown 
in Fig. 2(b), A C (L) behaves slightly convex upwards as a 
function of L^ 1 towards the bulk limit, A c (oo) < 10“ 2 , 
and the system is expected to be metallic. 

In conclusion, power law decay of the dimer-dimer cor- 
relation and the possible absence of a charge gap are con- 
firmed. At present, there is no evidence which denies the 
metallic property of the ground state at quarter-filling 
characterized by a Tomonaga-Luttinger liquid. 
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